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INTERACTION OF A TURBULENT-JET NOISE SOURCE WITH TRANSVERSE 
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MODES IN A RECTANGULAR DUCT 

by George P. Sued,* Kenneth J. Baumeister, and K. Uno Ingard^ 

Lewis Research Center 

SUMMARY 

A turbulent jet was used to excite transverse acoustic modes in a rectangular duct. 
The pressure spectrum showed asymmetric singularities (pressure spikes) at the reso- 
nant frequencies of the duct modes. This validates previously published theoretical 
results. These pressure spikes occurred over a range of jet velocities, orientations, 
and inlet turbulence levels. At the frequency of the spike, the measured transverse 
pressure shape matched the resonant mode shape. 

In the experiment, the pressure response of a duct to an enclosed air jet was 
measured at frequencies from 0 to 25 kilohertz and from 0 to 5 kilohertz. The jet 
noise in these experiments was produced by drawing air through an orifice into a 3. 81- 
by 10. 16 -centimeter (1^- by 4-in. ) rectangular duct. 


INTRODUCTION 

In the transmission of sound in a duct with or without flow, large amplitude varia- 
tions in pressure can be produced by the interaction of the sound with the duct walls. 

For example, the combustor noise signature in a turbojet engine (such as recently re- 
ported in refs. 1 and 2) can be modified by the ermine ducting. In those studies, the 
duct flow had some effect; nevertheless, the present study concentrates on the interac- 
tion of sound with duct walls in the absence of flow. Therefore, we used a turbulent jet 
to excite the transverse acoustic modes in a rectangular duct. Because the jet area 
was much smaller than the duct area, the duct flow was nearly zero and is thus ne- 
glected. The acoustic response to the jet near an eigenfrequency is identified, and the 

♦Massachusetts Institute of Technology, Cambridge, Massachusetts. 
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frequency range is determined for which the acoustic response in a duct differs qual- 
itatively from the response in free space. 

The qualitative acoustic response to turbulent jets in free space has lorg been known 
(ref. 3, pp. 123-125). In free space the sound pressure spectrum rises smoothly, ap- 
proximately as the frequency squared, to a peak at a Strouhal number (St = of 

about 0. 2 (ref. 4); it then falls off smoothly, approximately as the inverse of the fre- 
quency squared. In a duct, however, the pressure spectrum of a jet is more compli- 
cated than for the same jet in free space. 

Although the most popular model of the turbulent noise source of a free jet is a 
quadrupole (refs. 5 and 6), some success has been had in modeling the jet noise source 
as a distribution of simple (monopole) sources (ref. 7) . Also, a turbulent field in the 
presence of boimdaries has been shown to exhibit dipole source terms (ref. 6, ch. 3). 

If an enclosed jet can be modeled by a monopole source, the pressure spectrum theo- 
retically (refs. 8 and 9) will show asymmetric singularities (pressure spikes) at the 
duct-mode cutoff frequencies. 

Theoretically, the pressure spectrum of a fan noise source in a semi-infinite duct 
should show similar asymmetric spikes (ref. 6, p. 335). These spikes have been 
measured for a simulated fan noise source (ref. 10). Although they do not show the 
sharp spike at the resonant frequency predicted by theory, the measured spikes are 
asymmetric (ref. 10, fig. 8). In reference 10, however, detailed sound-pressure- 
level (SPL) traverses confirm that these peaks correspond to the excitation of higher 
order modes at the duct-mode resonant frequencies. Since both fan and jet noise 
sources can be modeled by simple acoustic sources, the simple source model should 
predict the response of a turbulent jet in a duct. 

The main purpose of this study was twofold: First, to confirm the existence of 
asymmetric resonance peaks in the acoustic spectrum of an enclosed turbulent jet, to 
show that the simple source model adequately predicts the resonant frequencies, and to 
identify the acoustic modes associated with resonance; second, to study the effects of 
jet velocity, orifice diameter and orientation, Strouhal number, and inlet turbulence. 
These experiments were performed in a new flow-duct facility; therefore, a secondary 
purpose is to document the facility and instrumentation. Also, some important aspects 
of pressure measurements peculiar to duct acoustics are discussed as an aid in inter- 
preting flow-duct data. 

The pressure response of a duct to an enclosed turbulent jet was measured at fre- 
quencies of 0 to 25 kilohertz. The jet noise was produced by drawing air through ori- 
fices of different sizes and orientations at different flow rates. 

SYMBOLS 

o 

Ah cross-sectional area of duct, m“ 
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radius of sphere about a simple source, m 

duct width, m 

speed of sound, m/sec 

orifice diameter, m 

duct height, m 

yTi 

wave number (see eq. (6)), m”^ 

axial propagation coefficient for mode number (m, n), m”^ 

source wave number of a propagating mode (see fig. 7 and eq. (18)) 

length of pipe, m 

Mach number, Vj^^/cq 

Mach number in pipe leadii^ to orifice 

transverse horizontal wave number 

number of propagating modes 

number of propagating modes in Ai/ 

transverse vertical wave number 

normal to wall 

acoustic pressure. Pa 

reference rms pressure. Pa 

rms pressure (see eq. (10)), Pa 

radiative resistance of source, g/sec 

Q 

source strength, m /sec 

Strouhal number, 

sound pressure level 

one period of oscillation, sec 

time, sec 

velocity, m/sec 

duct velocity, m/sec 

velocity of jet, m/sec 

transverse horizontal duct coordinate, m 
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transverse vertical duct coordinate, m 
axial duct coordinate, m 

Kronecker delta symbol (1 when i = j; 0 when i j) 
eigenvalue (transverse mode in duct), 1/m 
mean-square amplitude of cross -duct mode 
wavelength, Cq/v, m 
frequency, Hz 

cutoff frequency, Cq /2yWb?T(n/^, Hz 
Strouhal frequency, Hz 
frequency interval, Hz 

3 

atmospheric density, g/m 
wave function 

angular frequency, 2m’, rad/sec 


THEORY 

A sound wave propagating in a rectangular duct can be composed of a fundamental 
plane wave and various transverse waves that depend on the nature of the noise source. 
The fundamental plane wave, identified as the (0, 0) mode, travels in the axial direction 
and has a uniform pressure distribution across the duct cross section, A transverse 
wave, identified as the (m, n) mode, has a nonuniform pressure distribution across the 
duct; consequently, it reflects back and forth from the duct walls as it travels along 
the duct. 

The maximum number of modes in a sound wave propagating in a hard-wall duct 
depends on the frequency. For a given mode (i. e. , the (m, n) mode) to propagate, its 
frequency must be greater than the duct ’’cutoff” frequency. However, a noise source 
may not excite (at significant magnitude) the maximum number of modes that its fre- 
quency content would indicate. For example, a loudspeaker that fills the entire duct 
end might excite a large-amplitude plane wave but greatly reduced-amplitude higher 
order modes. In any case, as the frequency of the source is increased beyond the mode 
cutoff frequency, a new mode becomes available for the radiation of sound down the 
duct. 

The cutoff frequency for each mode and the pressure response of a simple (mono- 
pole) source in a duct are well known from classical acoustics (ref. 8, pp. 497-502). 
The classical theory relating to mode cutoff frequency and pressure response is briefly 
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reviewed here as a background for the experimental results. 


Cutoff Frequency 

Consider a rectangular duct of width b and height d. The wave equation describ- 
ing the propagation of an acoustic disturbance in a source-free region with negligible 
mean flow is 


V^P - — = 0 




( 1 ) 


where viscous effects have been neglected. 

Because of the hard walls the normal particle velocity must be zero, which is 
equivalent to requiring 


3P 

ant 


= 0 

walls 


( 2 ) 


where nt is the normal to the wall. The axes of the Cartesian coordinate system are 
chosen so that the walls are located at x = 0, b and y = 0, d, as shown in figure 1. 

For a sinusoidal time disturbance e^^^, the solutions to equation (1) are of the 
form 

P(x,y,z,t) ^ ^mn'^mn^^’ y)® 

m n 


where 
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and is an arbitrary constant. The functions are orthogonal to one another; 

that is, 


^3^ y)^pq(’‘> ?)*■ '*>' = WnqMnrn 

where is the cross-sectional area of the duct and the average value of 

The experiments were concerned with the excitation and propagation of the modes 
represented by the fimction as they travel down a duct. In particular, pressure 

traverses in the x-direction were made to determine if the ’^10’ ^^20 

were present. The pressure amplitude (spatial) distribution of the first three modes is 
shown in figure 2. The plus and minus signs on the duct cross sections in figure 2 in- 
dicate the spatial domain when the term takes on either a positive or a negative 

sign. (The term is defined in eq. (4). ) Recall from equation (3) that the pres- 

—“if 1+ 

sure P also varies harmonically, as e~ . The temporal variation of the mode 
is shown in figure 3. Although a microphone will respond to this instantaneous variation 
in pressure, a tsrpical measurement system records pressures in terms of the sound 
pressure level (SPL), defined as 


SPL = 10 


log 


10 


rms 


ref 


(9) 


where is the rms (root mean square) soiuid pressure defined as 

"■rms . O' 

-5 

where T is one period of oscillation. The reference rms pressure P^^^ is 2x10 
pascal. 

Assume, for example, that only one higher order mode is excited. Then substitut- 
ing the product of and the real component of e”^^^ into equations (9) and (10) 

yields the following SPL dependence on cross section: 
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SPL = Constant + 10 log 


( 11 ) 
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The pressure amplitude (spatial) and SPL variations of the first three modes are 
shown in figure 4. The minimum sound pressure levels in figure 4(b) correspond to the 
dashed partition lines in figure 2. These types of SPL distributions indicate the pres- 
ence of a particular mode. (The SPL is not, in general, equal to the intensity level in 
a duct if higher order modes are excited (ref. 9). ) 

We now replace the double index mn by the single index n by arranging the modes 
in order of increasii^ ^mn' equation (5), for the n " mode to propagate down 

the duct, must be real. For a given frequency v, is real only if (ref. 3, p. 504) 


” 2jt 


( 12 ) 


where is the cutoff frequency. That is, for any frequency, only a finite number of 
modes can propagate. To put it another way, each mode has a cutoff frequency: For 
frequencies below cutoff, the mode attenuates exponentially; for frequencies above cut- 
off, the mode propagates along the duct without attenuation. These ideas are illustrated 
in figure 5 for a modal pattern in a rectangular duct. The upper schematic shows decay, 
V < the lower schematic shows propagation, u > Dissipative mechanisms can 
cause some attenuation in the propagating wave. 


Impedance of a Simple Source 

The pressure near a simple source and the pressure radiated down the duct are 
directly proportional to the mechanical resistance at the source. However, the duct 
walls reflect sound back on the source and thus modify that resistance. The mechanical 
impedance of a point source inside a duct is discussed in the following paragraph. 

Consider the response of a medium in a duct of cross-sectional area to a 
simple (monopole) source of strength S and frequency v. The acoustic velocity is 
uniformly radial at a spherical surface of radius a (a « X) that just encloses the gen- 
erating mechanism. The mechanical resistance of such a source in a duct can be shown 
to be (ref. 8, p. 502, eq. (9. 2. 12) 


R = 


Pq^O 


Q 2 4 
on a 
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(13) 
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where the source is at the cross-sectional coordinate point (xq, Yq). The summation in 
equation (13) includes only the propagatii^ modes, the last of which mode is labeled N. 

To gain some insight into equation (13), we rewrote it for three ranges of the fre- 
quency V (recall that represents the cutoff frequency): 


For 0 <. V < 


R = 


2 4 

PqCqStt a 


(14) 


For V < v^. 


R = 


2 4 


^ , k ^i(vV 

k. A, 


(15) 


For 1^2 — ^ < ^3* 


R = 


2 4 

PqCqStt a 


kf Aj k2 A 2 


(16) 


and so forth, as v increases. 

A discontinuity (pole) appears in R as each cutoff frequency is crossed, since 



for V - u^. As V increases beyond the mechanical, or equivalently the radiative, 
resistance R drops as k^^. When u increases above a new eigenfrequency, another 
pole is included in the expression for the resistance. 

Figure 6 illustrates these resistance characteristics for. an example taken from 
reference 8. As shown in figure 6, when the frequency of a simple source in a duct 
just equals the cutoff frequency, the theoretical radiation resistance increases step- 
wise and then decreases slowly as the frequency goes beyond cutoff. The number by 
each peak is the mode number (defined in eq. (4)). The dashed curve is the radiative 
resistance of the same source in free space. 
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Pressure Response of a Simple Source 


As mentioned, the pressure on the source and the pressure radiated to the far field 
are directly proportional to the resistance R. However, real fluid effects - such as 
energy dissipation by viscous friction, and thermal conduction - limit the pressure re- 
sponse at a cutoff frequency to a finite value and increase the bandwidth. Nonetheless, 
discrete jumps should occur at the cutoff frequency. 

As the frequency of the source increases (eq. (13) and fig. 6), mechanical resist- 
ance increases in discrete jumps, each jump corresponding to an additional propagating 
mode. As cu becomes large, the ninnber of propagating modes added per unit change 
in 00 becomes larger, the interval between modes becomes shorter, and the peaks 
become narrower. At the high-frequency limit (ref. 8, p. 503), 

limit (R) =p„c„(47ra^)k^a^ (17) 

00 ^°o ^ ^ 


The average resistance is the radiative resistance of a source in free space. Conse- 
quently, at high frequencies the pressure response of a source in a duct approaches the 
response of the same source in free space. 

For example, consider a rectangular duct - the geometry used in these experiments. 
The high-frequency limit of the radiative resistance is based on the number of modes 
available for sound propagation. Therefore, counting these modes gives some insight 
into their distribution with frequency. 

The central effect is the phenomenon of mode cutoff. Cutoff can be discussed in 

terms of the characteristic mode wave numbers . Recall that the wave number 

m, n 

is the angular frequency co divided by the speed of sound Cq. A mode propagates (does 
not exponentially decay) only if the source wave number kg = oo/cq is greater than or 
equal to ^ (cutoff value of k) as given in equation (7). 

The problem is illustrated graphically by figure 7. The mode wave numbers are 
indicated by dots on a wave number graph (fig. 7(a)) and by a bar on the wave number - 
frequency plot (fig. 7(b)). For a mode to propagate, it must be included in the pie- 
shaped region described by radius kg (fig. 7(a)). The total ntunber of modes in this 
region is given by the number in the center of the pie -shaped region plus the number 
along each axis (ref. 8, p. 586), which is approximately 
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It is also useful to note how many new modes are available per unit increase in 
wave number Akg. Graphically, the number of modes added at a given frequency is 
proportional to the area of the quarter annulus surrounding the pie-shaped region. Ob- 
viously, this area is directly proportional to the source wave number kg. From equa- 
tion (18), the number of new modes that can propagate N per unit increase in kQ is 

AN N(kQ + Akg) - N(kQ> ffkQ ^ j k^ ^2 ^ 

= — ~ 4 1 = 

Ako AkQ 2 (-~\ (-\ ^ ^0 (--\ 

\b d/ \b/ Vdj \b d/ (^b d/ 

At small wave numbers (frequencies), AN/Akp is small; thus the frequency inter- 
val between the cutoff mode and the next propagating mode is large. At large wave 
numbers, AN/Akg is large; thus the frequency interval between the cutoff mode and 
the next propagating mode is small. Equation (19) shows that AN/Akg is directly pro- 
portional to frequency. 

Experimental evidence of pressure spikes in a duct can be found only in the low- 
frequency range. Here the separation between modes is great enough to exhibit the 
asymmetric shape of the pressure spike. In the high-frequency range, neighboring 
modes overlap and obscure the shape of the pressure spikes. Pressure spikes are not 
observed at high frequencies since the bandwidth of a given pressure spike is much 
wider than the frequency interval between modes. The bandwidth of a spike (fig. 6) 
may be widened by energy dissipation. Thus, spike overlap is a consequence of the 
duct’s acoustical properties. Increasing the frequency resolution of the detecting sys- 
tem therefore may not reveal any pressure spikes. 


Theoretical Verification 

Various theories have been postulated for the mechanisms by which jets generate 
sound: for example, the quadrupole, dipole, and monopole theories. In a recent pres- 
entation (abstracted in ref. 11), Powell pointed out that the various theories have lead 
to interesting and apparently divergent views of the noise-generating source. However, 
regardless of the source, the present experiment should have revealed asymmetric 
pressure spikes at the ’’low" cutoff frequencies and a blending of the pressure into the 
free-space pressure response at the "high" frequencies. The reason is that pressure 
spikes depend primarily on the propagation characteristics of the medium and only 
secondarily on the noise source. For example, a dipole can be considered as two 
monopoles, 180° out of phase, separated by a small distance. Each monopole radiates 
as if the other was not there, and the radiative field is the sum of the two - asymmetric 
spikes and all. 
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APPARATUS and PROCEDURE 


To verify the theory, we first placed a simple, broadband, noise source capable 
of exciting all the modes inside an infinitely long, hard-wall duct and then monitored 
the pressure spectrum for the characteristic asymmetric pressure spikes at the duct 
cutoff frequencies. The flow-duct system shown in figure 8 was adapted to simulate 
this idealized infinite-duct situation. The following paragraphs described this flow-duct 
system; first, the general operation; then, the key system components, in detail. 

As shown in figure 8, a flow duct was connected by a 15-centimeter- (6-in. -) 
diameter vertical pipe into the NASA Lewis 1. 2-meter- (48-in. -) diameter altitude ex- 
haust header. The exhaust system has a nominal 68 940-pascal (20-in. -Hg) pressure 
at this point. The system was placed on line by first activating the 15-centimeter- 
(6-in. -) diameter remote-electric-motor valve directly beneath the 1. 2-meter- 
(48-in. -) diameter header line and then opening the pneumatically adjustable, 
15-centimeter- (6-in. -) diameter butterfly valve at the upper end of the vertical pipe. 

Sound is generated by pulling room air through a small inlet orifice in the test 
section. The advantage of this suction system is that any turbulence in the flow is in- 
duced solely by the inlet orifice. In constrast, a compressed-air system would require 
a large settling tank to insure uniform flow conditions. And a settling tank would re- 
quire acoustic treatment to prevent the excitation of resonant modes (ref. 12). 

In the suction system the air jet passing through the orifice into the duct is suffi- 
cient to generate a strong (80 to 130 dB) noise source with the characteristic jet spec- 
triim. However, the small-orifice jet diffuses into the large duct in such a way that the 
duct Mach number effects are negligible. The flow rate was set by adjusting the 
15-centimeter- (6-in. -) diameter butterfly valve and placing porous resistance in the 
duct just upstream of the test section. 

The duct was made to appear acoustically infinite in length by placing absorbing 
material both upstream and downstream of the test section. The pressure response of 
the duct was measured by flush, wall-mounted microphones. 


Flow Duct 

The rectangular duct has a 3. 81-centimeter by 10. 16-centimeter (1^-in. by 4-in. ) 
cross section (fig. 9). An exponential horn was attached to the upstream end of the test 
section. The horn can be removed if the experiment dictates. For these experiments 
the horn was filled with acoustic foam to approximate a PqCq termination, which simu- 
lates an infinitely long duct. 

The 1. 4-meter- (4-ft, 8-in. -) long section just downstream of the horn is the 
actual test section. It consists of 16 detachable plates (8 on top and 8 on the bottom). 
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These plates are used to mount microphones, orifices, or pressure probes as the ex- 
periment dictates. 

The next 1. 8 meters (6 ft) of the duct contain an acoustic muffler made from two 
strips of acoustic foam (3. 8 cm by 3. 2 cm by 183 cm; 1-| in. by 1-| in. by 6 ft). These 
strips are mounted on a wooden support frame and inserted into the duct. The muffler 
serves two purposes: First, to attenuate any noise generated by the suction apparatus 
downstream of the test section; and second, to minimize the axial duct modes caused by 
reflection from the duct ends and by discontinuities of the duct. The muffler's trans- 
mission loss exceeds 40 decibels for frequencies above 1 kilohertz and 20 decibels for 
frequencies above 250 hertz. 

Fixed to the end of the rectangular duct is a circular pipe that leads to the suction 
apparatus. As mentioned, the flow is controlled by two valves: A cutoff valve connects 
the apparatus to the main exhaust facility; a pneumatically operated butterfly valve 
modulates the mass flow rate. 

The existing valves were designed to handle mean duct flows with Mach ninnbers of 
about 0. 5. In these experiments the mass flow rate was set by the desired flow through 
the orifice, which corresponds to duct Mach numbers of about 0. 01. Unfortunately, 
leaks in the butterfly valves (tolerable in high-mass-flow-rate experiments) produced 
excessive mass flow rates. To overcome this problem, we installed a flow resistor, a 
wire screen in front of which we placed blocks of acoustic foam, between the muffler 
and the test section. Since varyir^ the number of foam blocks varied the flow resist- 
ance and lowered the pressure drop across the orifice, we achieved finer control of the 
jet velocity. 


Noise Source 

Drawing air through a small orifice generates a turbulent jet, which will excite all 
the duct modes. Each duct mode will decay or be amplified depending on how it reacts 
with the duct. Therefore, the turbulent jet is ideal for studying the interaction of a 
noise source with the duct walls. In contrast, for example, an electrically driven 
speaker acts like an oscillating piston whose noise source is made up predominately of 
a plane pressure wave that will not react with the duct walls. Circular orifices were 
used to force the air to enter the duct as a round jet. Three orifices 0. 64, 1. 27, and 
2. 54 centimeters (1/4, 1/2, and 1 in. ) in diameter were machined into a 0. 32- 
centimeter- (1/8-in. -) thick plate (fig. 10(a)). The holes were countersunk to present 
a sharp edge to the incomir^ flow. This precaution was necessary because cylindrical 
holes with thickness -to-diameter ratios between 1/2 and 2 exhibit screech tones, which 
we want to avoid. In addition, a 1. 27 -centimeter- (1/2-in. -) diameter tube 1. 5 meters 
(5 ft) long was fitted to a test plate (fig. 10(b)) to produce a different type of jet noise. 
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Both the orifice plates and the long tube were mounted on top of the duct so that the jet 
would enter perpendicular to the duct axis, as shown in figure 11. The thin-orifice 
jet is initially laminar; the long-tube jet is initially turbulent. 

To allow the air to enter parallel to the duct axis, we replaced the exponential horn 
with a plate, shown in figure 12, that contained a 1. 27-centimeter- (1/2-in. -) diameter 
orifice. Later, results from this configuration are compared with results from the 
configuration shown in figure 11. 

In these experiments, the orifices were mounted in the test section so that the 
countersunk side faced the inside of the duct and thus presented a smooth surface to the 
incoming flow. Consequently, there were no mechanical variations at the entrance to 
influence sound generation. 


Flow Measurement 

The system orifice plate shown in figure 8 has a 9. 78 -centimeter- (3. 85-in. -) 
diameter hole. This plate was sized to measure flows with Mach numbers of about 0. 5; 
consequently, it could not be used to measure the small flows encountered in these ex- 
periments. Instead, the pressure drop measured across the orifices through which the 
flow is drawn into the duct was used to estimate the jet velocities. 

The orifice velocities were estimated from the atmospheric pressure and the static 
pressure just upstream of the exponential horn. Standard orifice equations (ref. 13, 
section 4-61) were used to calculate the orifice jet velocities for a discharge coefficient 
of 0. 61 and an adiabatic process. These calculated velocities, which are presented 
later with the acoustic results, are only rough estimates of velocity, not exact values, 
since the flow streamlines and instrumentation do not match the standard orifice con- 
figuration. 

As an additional cross-check, a hot-wire anemometer was used to measure the 
mass flow rate for both end- and side-oriented orifices (fig. 13). The square of the 
probe voltage is approximately proportional to the square root of the mass flow rate. 
The mass flow rate is predicted from the pressure drop across the orifice in the man- 
ner described. Each value of mass flow rate corresponds to a measured probe voltage. 
Since the flow rates for the side- and end-oriented orifices coincided, the flow rate ap- 
parently is not sensitive to orifice orientation. Consequently, the standard orifice 
equations probably give reasonable estimates of the flow rate. 


Microphone Procedure 

In an e3q)erimental run, the pressure in the duct was reduced to below atmospheric. 
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Air then flowed through the orifice, hit the opposite wall (fig. 11), and flowed down the 
duct into the altitude header (fig. 8). The turbulent fluctuations generated by the ori- 
fice jet acted as an acoustic source. 

These acoustic fluctuations were detected by a microphone at a certain minimum 
distance from the orifice. Two effects determine this minimum distance: 

(1) Near the jet, in the mixing region, the turbulent fluctuations are stroi^ enough 
to completely mask the acoustic fluctuations; therefore, the microphones must be out- 
side this region. In these experiments, the microphone was upstream of the orifice 
near the exponential horn (fig. 9), where the mean flow is zero. Here, the convected 
turbulence is negligible. Data are presented later on to assess the effect of micro- 
phone placement (see section Effects of Microphone Location). 

(2) The second consideration involves the detection of higher order acoustic modes. 
Although the frequency of a disturbance is below the mode cutoff frequency, the trans- 
verse pressure waves associated with the mode will exist near the source. These 
pressure waves, however, decay exponentially with distance. If the lowest mode is to 
be attenuated by a factor 1/e, the microphone should be placed 2b from the source, 

b being the width of the duct, where it will detect only propagating modes. For 
frequencies just below cutoff, the cutoff mode may be detected depending on how close 
the frequency is to the cutoff frequency. 

In the experiments, two 0. 64-centimeter- (1/4-in. -) diameter condenser micro- 
phones were used to determine the acoustic field. The response of these microphones 
was flat to 25 kilohertz. A standard acoustic source was used to check the system 
amplification. This source generated a 114-decibel signal in 5-octave steps from 
0. 125 to 2. 0 kilohertz. The microphone was placed in the source and the resulting 
decibel level of the output was noted. In this way an absolute scale was established for 
the SPL. One microphone was mounted with its membrane flush with the surface of a 
test plate (protective shields were removed). A typical microphone position is shown in 
figure 14. In general, only one microphone was used for most tests, except in one ex- 
periment comparing up- and downstream radiation patterns. 

One purpose of the present experiment was to detect the excitation of propagating 
higher order duct modes. Higher order acoustic modes are characterized by their 
pressure variations across the duct cross section. For it to detect these modes, the 
microphone must be at a pressure antinode, which in a rectangular duct is in a corner. 
Two plates were made with the microphones mounted in a corner. A third plate was 
machined to allow the transverse position of the microphone to vary (fig. 14). This 
plate had a double row of 0. 64-centimeter- (1/4-in. -) diameter holes with 1. 27- 
centimeter (1/2-in. ) separation between the centers. The second row was staggered 
behind the first to vary the microphone location by 0. 64 centimeter (1/4 in. ). The 
holes were stopped with aluminum plugs when the microphone was not in place. 
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Instriimentation 


The measuring apparatus used in this study is shown in the block diagram of fig- 
ure 15. The acoustic results were obtained from the x-y recorder, which recorded 
the output of the spectrum analyzer. Two frequency ranges were explored: 0 to 5. 12 
kilohertz and 0 to 25. 6 kilohertz. Each range consisted of 256 frequency bands to pro- 
vide bandwidths of 20 and 100 hertz, respectively. The temperature in the duct was 
assmned to be ambient. 

The 0- to 5-kilohertz frequency range spans the first four higher order modes in 
the duct. In this range the frequency separation between modes is sufficient to reveal 
the excitation of the individual modes. In the 0- to 25-kilohertz range, however, the 
modal density AN/AkQ, defined by equation (19), increases with frequency. In this 
range, the spectrum analyzer will not be able to discern discrete pressure spikes at 
the higher cross-mode eigenfrequencies. 


Summary of General Procedure 

In each experiment, an orifice was mounted in the test section, and then the flow- 
duct system valves were opened to reduce the pressure inside the duct. The flow rate 
through the orifice was set by inserting flow resistance in the duct as required and 
adjusting the butterfly valve. A microphone was placed in the far field of the source, 
and the instrumentation was activated to record the noise signal. The noise signal was 
Fourier analyzed, and separate averages of 256 channels were computed by the spectrum 
analyzer and plotted on the x-y recorder. 


DISCUSSION OF RESULTS 

In this section, the experimental results are discussed in the following categories: 
basic effects and the effects of jet velocity, screech, Strouhal frequencies, inlet tur- 
bulence, duct resonances, axial orientation of orifice, microphone location, and mode 
shape . 


Basic Effects 

In the basic experiment, a 1. 27-centimeter- (1/2-in. -) thick orifice plate was 
mounted in the test section as shown in figure 9, and air was drawn through it into the 
duct. The turbulent fluctuations thus generated act as an acoustic source. The result- 


ii^ acoustic pressure oscillations travel through the duct and are detected by the micro- 
phone system. The noise spectra associated with these oscillations are shown in fig- 
ure 16. 

Recall from earlier theoretical discussions that the influence of modes on the 
acoustic spectrum is most evident at low frequencies, when the separation between 
modes is large and the nvunber of propagating modes small (fig. 7). At high frequencies 
the effect of an individual mode is averaged out; hence the spectrum tends toward that 
of a jet in free space. In the 0- to 5-kilohertz range, the SPL increases abruptly near 
the mode cutoff frequency (fig. 16(a)). The effect is quite pronounced in the first three 
modes. In the 0- to 25-kilohertz range, the acoustic spectrum changes character 
(fig. 16(b)). Above 7 kilohertz, pressure spikes are no longer clearly evident at the 
eigenfrequencies. Pressure spikes of neighboring modes overlap and minimize the in- 
fluence of any single mode. The modal density eventually becomes too great for the 
256-channel spectrum analyzer to separate the individual pressure spikes; consequently, 
the spectrum represents an average of many modes. To this extent, the SPL spec- 
trum above 7 kilohertz (fig. 16(b)) is roughly comparable to the SPL spectrum of a 
jet exhausted into free space, which confirms the theory. 

Thus, the basic experiment had the anticipated results. 


Effects of Jet Velocity 

To determine the effect of jet velocity on the acoustic spectrum, we chose several 
pressure drops across the orifice. Representative acoustic spectra over a range of 
jet velocities are presented in figure 17 for the 0. 64-, 1. 27-, and 2. 54-centimeter- 
(1/4-, 1/2-, and 1-in. -) diameter orifices mounted on the side of the plate as shown in 
the sketch in figure 17(a). 

The spectra for the 0- to 5-kilohertz range (fig. 17) reveal pressure spikes at each 
mode eigenfrequency (see also fig. 16). These asymmetric pressure spikes persist 
throughout the range of orifice diameters and Mach numbers tested. However, the 
overall acoustic level rises as the Mach number is increased. Moreover, the spectral 
intensity shifts to the higher modes as the mean velocity is increased, which is shown 
by the magnitudes of the first and fourth modes in figures 17(a) and (c) for increasing 
Mach number. Figure 17(b), however, does not follow this trend, particularly at 
Mach 0. 44. The SPL of these fluctuations increases with increasing velocity. Thus, 
for a source of this type, there is generally an increase in level with velocity as well 
as a shift in level to higher frequencies. 

The spectra in the 0- to 25-kilohertz range (fig. 18) are relatively smooth, like a 
free-jet spectrum, above 5 kilohertz. In figures 18(a) and (b), dashed lines repre- 
senting the average behavior of the spectra are drawn through the highest velocity 
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curves. In figure 18, the peak Strouhal number is about 0. 2, which is like that of a 
free jet (ref. 4). 


Effects of Screech 

Screech peaks occur at the higher Mach numbers in figure 18. These peaks are not 
caused by the action of the duct walls on the jet noise source. If they were, the curve 
would resemble that for low-frequency oscillation - a discontinuous jump at a cutoff 
frequency and a slow decay in level above the cutoff frequency. Since the screech 
peaks are symmetric about their maximum, there is probably a different cause. 

Screech is generally considered a nonlinear process in which the acoustic field can re- 
act back on the jet. Under this condition, intense pure-tone oscillations that are caused 
by a jet impinging on a plate can be detected (ref. 14, ch. 4). 

Similar symmetric noise peaks also occur, for example, when air flows over 
resonator cavities (refs. 15 and 16). The shedding of vortices by the flow over the 
resonator cavity holes coupled with the cavity resonance or even the duct mode reso- 
nance (ref. 15) produces these symmetric noise peaks. For screech and cavity- 
induced noise, the simple source model does not adequately describe the feedback 
mechanism. 


Effect of Equal Peak Strouhal Frequencies 

The average shape of the spectrum is dominated by the location of the peak fre- 
quency. The distribution in acoustic intensity with frequency should be similar for dif- 
ferent orifices with approximately equal peak Strouhal frequencies. However, the ab- 
solute value of the intensity will be higher for the larger diameter orifice with the 
higher jet velocity, as shown in figure 19 for three orifice diameters and three Strouhal 
frequencies. 

To look in more detail at the rise in overall acoustic level with increasing velocity 

o 2 

and diameter, we assumed that the radiative power in the duct is proportional to v D 
(ref. 6, p. 93), as for free jets. If the spectrum shapes are identical, the power per 
unit frequency should also increase. The acoustic level ratio between orifices of differ- 
ent diameters but similar Strouhal frequencies is 


AS PL = 10 log 





or 


Vj D. 

ASPL = 80 Ic^ — + 20 Ic^ — 


'0 


D, 


0 


For equal Strouhal numbers, (j^Dq/vq) equals (i/Dj/v^); therefore, if the orifice 
diameter is doubled, as in figure 19, the velocity must double to keep the peak Strouhal 
number constant; thus 


D, 


Dr 



2 


and 


ASPL = 80 log 2 + 20 log 2 = 30 dB 

The difference in level of both spectra in figure 19 is approximately 30 decibels, as 
predicted. 

Although this result holds for high-frequency spectra, where the duct response is 
similar to the free-space response, it is only approximate for low-frequency spectra. 
One noteworthy exception is the behavior of the spectra below the first cutoff frequency. 
This behavior is most pronounced in figure 19(c) below 2. 0 kilohertz, where the spec- 
trum for the 0.64-centimeter- (1/4-in. -) diameter orifice with the lower Mach number 
decays from = 0 to the SPL pressure spike at the first cutoff frequency. However, 
the spectrum for the 1. 27-centimeter- (1/2-in. -) diameter orifices with the higher 
Mach number does not decay. (Note that the low-frequency spectra of fig. 17 also 
decay with frequency at low Mach numbers. ) 

At low Mach numbers the behavior below the first mode may be dominated by the 
sound resulting from the complex recirculation patterns along the wall (fig. 11) rather 
than by the orifice jet noise, as shown by the schematic in figure 20. That is, as the 
level of the jet spectrum is decreased, the recirculation noise source could be un- 
covered, which could account for the dip in the SPL spectrum before first-mode cut-on. 
The spectrum shape below the first cutoff will, of course, depend on the relative 
strengths of the various duct noise sources. 


Effect of Inlet Turbulence 

If the pressure spikes in the 0- to 5-kilohertz frequency range are due solely to 
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the duct's response to a turbulent jet source, changes in inlet turbulence should not 
obscure this effect. To test this thought, we altered the turbulent nature of the inlet 
jet flow. The flow entering the duct from a square-edge orifice is laminar; it becomes 
turbulent further downstream. The flow entering the duct from a long tube, on the 
other hand, is fully turbulent. If the duct's response depends only on the turbulent jet 
sources, the change in the jet flow from laminar to turbulent near the tube outlet should 
not destroy the effect. These considerations were checked by allowing the flow to enter 
from a 1. 5-meter- (5-ft-) long tube of 1. 27-centimeter (l/2-in. ) inside diameter mounted 
on the side of the duct (fig. 12). 

The acoustic spectra in the 0- to 5 -kilohertz range (fig. 21(a)) still contain discrete 
pressure spikes at the cross-mode frequencies. Thus, the initial turbulence level of 
the jet core probably has a very small effect on the turbulent shearing region. The 
acoustic spectra in the 0- to 25-kilohertz range (fig. 21(b)) again resemble the free- 
space spectra at high frequencies. A peak in level occurs at approximately the peak 
Strouhal number of 0. 2, which is consistent with previous behavior. Screech peaks are 
absent from these curves; however, tlie velocity range explored is not sufficient to re- 
solve the influence of inlet turbulence on screech. Screech has been observed in turbu- 
lent jets (ref. 14, section 4. 3. 5). 


Effect of Duct Resonances 

Periodic oscillations in amplitude with frequency resulting from resonance in the 
long tube appear in figure 21(a). Acoustically, the inlet tube acts like open organ pipes, 
for which the resonant frequencies are given by 


( 2 ») 

where the wavelength of the sound must be larger than the tube diameter, and the in- 
ductive end correction is neglected. The s3mibol represents the average flow Mach 
munber in the tube. The difference in frequency between two adjacent resonances of the 
pipe is 




n+1 


n 




( 21 ) 


For a 1. 5-meter- (5-ft-) long pipe with Mach 0. 2 flow, At> is 108 hertz, which is close 
to the measured spacing between peaks of 112 hertz. Therefore, the ripples in ampli- 
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tude with frequency in figure 21(a) represent the excitation of the axial pipe modes that 
are superimposed on the duct’s acoustic response. The peaks in the ripples represent 
the axial resonance condition in the pipe. 

Similarly, axial (z-direction) acoustic modes can be excited by sound-wave reflec- 
tions at the ends of the flow duct. In the study of reference 17, symmetric pressure 
oscillations were induced in a fan inlet by varying the inlet length to induce axial reso- 
nance. In the present study, axial resonances were minimized by the entrance muffler 
and the foam -filled exponential horn (fig. 9). 


Effect of Axial Orientation of Orifice 

The position of the orifice was altered so that the jet entered the duct parallel to 
the duct axis: The exponential horn was removed from the upstream end and a 1. 27- 
centimeter- (1/2 -in. -) thick orifice plate was inserted in its place (fig. 13). The duct 
then resembled a semi-infinite duct with a hard-wall entrance, rather than the infinite 
duct developed in the theory. 

Acoustic spectra in the 0- to 5-kilohertz range as a function of Mach number 
(fig. 22) again show the basic effects of discrete level changes at the cross-mode eigen- 
frequencies. The noise levels for side- and axially oriented holes are about the same 
(see also fig. 17(b)). However, the considerably higher low-frequency noise for the 
axial jet tends to obscure the first peak. Because of this and because the side-oriented 
orifice gave results more nearly matching theory, we chose to use side-oriented ori- 
fices. 


Effects of Microphone Location 

Except for experiments with the axially oriented orifice, all experiments were per- 
formed with the microphone located upstream. In the upstream region, there is no 
convected turbulence to mask the acoustic fluctuations. In addition, there is no change 
in duct response with jet velocity, because there is no mean flow upstream of the ori- 
fice. However, since the mean flow downstream of the orifice is low, there should be 
no pronounced differences in the upstream and downstream radiation. 

To test the effects of microphone location, we used two microphones: one upstream 
and one downstream of a side-oriented orifice. The acoustic spectra from the two mi- 
crophones were then plotted on the same graph for jet Mach numbers of 1. 00 and 0. 25 
(fig. 23). At Mach 0. 25 the responses of the upstream and downstream microphones 
differed by only 0. 8 decibel. At Mach 1. 00 the difference was greater: The downstream 
spectrum varied periodically in level with frequency. And, though the downstream 
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spectrum was, on the average, identical to the upstream spectrum, it somewhat ob- 
scured the duct's response to a turbulent source. To avoid this, we placed the micro- 
phone upstream in all other experiments. 

Another test was run while a second microphone was traversed axially along the 
centerline of the duct. As shown in figure 24, the position of the traversing microphone 
changes the SPL level detected by the stationary wall microphone and thus changes the 
standing wave pattern in the duct. Reflections from the traversing microphone domi- 
nate the duct response, most likely because of the relatively small duct cross-sectional 
area. Thus, insertii^ any probe into the flow duct of this facility may significantly 
alter the duct acoustics. 


Effect of Mode Shape 

These experiments are intended to show conclusively that the duct response is due 
to the excitation of the higher order modes and to emphasize that the microphone must 
be placed at a pressure antinode. Therefore, the pressure level in the x-direction 
(fig. 1) was measured at three cross-mode eigenfrequencies: 


Frequency, 

Number of 

Propagating 

kHz 

propagating 

modes 


modes 


1 

1 

(0,0) 

1. 8 

2 

(0,0), (1,0) 

3. 5 

3 

(0,0), (1,0), (2,0) 


As the frequency is increased above the mode cutoff frequency, the previously excited 
modes still propagate. However, slightly above the mode cutoff frequency, the duct 
response is dominated by transverse (x-direction in this case) variations of that reso- 
nant mode. The reason is the pole in the duct response (eq. (10)) at the cutoff fre- 
quency. 

The measured SPL variations of the (0, 0), (1, 0), and (2, 0) modes shown in fig- 
ure 25 correspond to the theoretical variations shown in figure 4. At 1 kilohertz, only 
the plane wave propagates and pressure does not vary across the duct cross section. 

At 1. 8 kilohertz, the first asymmetric mode dominates the response - note the node at 
the center. At 3. 5 kilohertz, the first symmetric mode dominates the response - note 
the antinode at the center. 

The pressure fields in the various modes are additive. If a microphone is placed at 
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a pressure node for a given mode, the mode will no longer be detected. For the first 
asymmetric mode, we compared the acoustic spectrum at centerline to that at the 
corner (fig. 26). Because the center microphone position is a node for this mode, 
the peak was not detected. 


CONCLUSIONS 

The following conclusions have been reached about a turbulent noise source in a 
rectangular duct: 

1. In the frequency range 0 to 5 kilohertz, which includes only the first few duct 
modes, the duct walls have a strong influence on the acoustic spectra. The discrete 
nature of the first few higher order modes manifests itself in an as 3 nnmetric singularity 
(pressure spike) near the mode cutoff frequency. 

2. In the frequency range 0 to 25 kilohertz, where the number of propagating modes 
is large and the intervals between modes are small, the influence of the duct walls is 
minimal and the acoustic spectrum is nearly as smooth as if the jet were in free space. 

3. The asymmetric pressure spikes in a duct occur over a range of jet velocities, 
jet orientations, and inlet turbulence levels. 

4. The qualitative structure of the frequency spectra can be interpreted through a 
linear analysis of the duct's acoustical response. 

8 2 

5. For enclosed jets, the intensity level appears to be proportional to v D , where 
v is velocity and D is orifice diameter. That is, the relationships for free jets apply 
equally well to enclosed jets, for the limited data taken. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 24, 1978, 

505-03. 
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Figure 1. - Duct coordinate system. 




(cl Mode 0). 

Figure 2. - Mode concept - spatial variation in pressure. 


Amplitude 



(a) Elapsed time, 0. (b) Elapsed time. T/8. (c) Elapsed time, T/4. 



(d) Elapsed time, 3T/8. (e) Elapsed time. T/2. (f) Elapsed time, 5T/2. 



(gl Elapsed time. 3T/4. (h) Elapsed time, 7T/8. (I) Elapsed time, T (one 

period of oscillation). 

Figure 3. - Time dependency of spatial pressure variation for (1. 0) mode. 



(a-3) Mode (2, 0). (b-3) Mode ( 2 , 0). 

(a) Spatial pressure variation. (b) SPL pressure variation. 

Figure 4. - Spatial pressure and sound pressure level (SPL) variations for 
modes (0,0), (1,0), and (2,0). 
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Exponential amplitude decay with x 


la) Sound frequency below cutoff frequency, w < u 


Propagation 
down duct 


(b) Sound frequency above cutoff frequency, v > 

Figure 5. - Illustration of decaying and propagating modes, where u is mode 
oscillation frequency and is cutoff frequency (predicted by theory). 


-Radiation resistance 
to same source in 
free space 


Figure 6. - Dimensionless radiative resistance R of a simple source of radius a at point (XQ,yQ) in a 
rectangular duct of sides b and d, where Xq = 1/4 b; yg = 1/10 d; and d = 4/5 b. (From ref. 8. ) 
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(a) Wave number graph for duct whose width is twice its height. 
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(b) Mode diagram for rectangular duct 3.81cm (1.5 in.) high by 10.16 cm (4 in.) wide. 
Rgure 7. - Graphic illustration of mode cutoff. 
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Figure 9. - Flow-duct test section. 
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Orifice 
diameter, 
cm (in. ) 

0.64(1/4) 2.54(1) 1.27(1/2) 



(a) 0. 32-Centimeter-thick plates. 



(p) Long tube. 

Figure 10. - Orifice plates. 
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Figure 11. - Flow pattern of jet into duct. 



Figure 12. - Flow-duct apparatus. 
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Figure 15. - Instrumentation. 



■g (a) Low-frequency range. 

O 



0 5 10 15 20 25 


Frequency, v, kHz 
(b) High-frequency range. 

Figure 16. - Low- and high-frequency spectra in a rectan- 
gular duct. Orifice diameter, 1.27 cm (1/2 in. h orifice 
orientation, top of duct; Mach number, 0.44. 
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Sound pressure level, dB (referenced to 2(10"^ Pa) 


jgj o Micropnone 


Mach 

number 



'>-Screech 

peaks 



Frequency, v. kHz 


(a) Orifice diameter, D, 0.64 cm (l/4in. ). 


(b) Orifice diameter, D, 1. 27 cm (1/2 in. ). 
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(c) Orifice diameter, D, 2.54 cm (1 in.). 

Figure 18. - Effect of orifice jet velocity (Mach number) and orifice diameter on acoustic spectra in the high-frequency range 
(Oto 25 kHz). 
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Frequency, v, kHz 


(a) Strouhal frequency, Uj, -1.9 kHz. (b) Slrouhal frequency, Uj, -2.5 kHz. 



Frequency, v, kHz 


(c) Strouhal frequency, Vj. -5.7 kHz. 

Figure 19. - Sound pressure levels for orifices with approximately equal peak Strouhal frequencies (Uj = 0. 2 v/D). 




Frequency — ^ 

Figure 20. - Schematic representation of noise sources in flow duct. 



Frequency, v, kHz 

(a) Low-frequency spectra. (b) High-frequency spectra. 

Figure 21. - Effect of velocity of flow entering orifice from long tube on acoustic spectra. Tube length, 1. 52 m (60 in. h tube 
inside diameter, 1.27 cm (l/2in. ). 
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Figure 2Z - Low-freguency spectra observed when jet and 
duct axis are parallel. Orifice diameter, D, 1.27 cm 
(1/2 in.). 
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Figure 23. - Upstream- and downstream-radiated acoustic 
spectra. 



Relative sound pressure level, dB 



Figure 24. - Minimuin-to-peak sound pressure level variations as 
function of microphone location. 



Microphone iocation, x/L 


Figure 25. - Transverse variation of radiated sound ievels at frequencies seiected to show mode shape. 
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